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0\ . ABSTRACT 
^ ' 

We investigate a problem of the necessary and sufficient conditions for ap- 
pearance of the 1/f fluctuations in the simple systems affected by the external 
random perturbations, i.e. the power spectral density of the flux of particles moving 
in some contours and perturbed by the external forces. In some cases we observe 
J> ! the 1/f 5 (8 ~ 1) behavior but only in same range of frequencies and parameters of 
^ | the systems. 
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INTRODUCTION 

O ■ 
00 ■ 

Q\ ■ Many processes of the greatest diversity exhibit temporal evolution charac- 

(3JT)' terized by a power spectral density S(f) which behaves like 1/f 5 (5 ~ 1) at low 
frequencies. Despite considerable efforts, the generally accepted theory of the 1/f 
noise is still lacking. Possibly, there are several physical mechanisms giving rise to 
1/f noise. 

Usually a study of 1/f noise in materials starts from the approximate em- 
pirical relation given in 1969 by Hooge which for the current / may be written as 

[i] 

Here I is the squared average current, a is a dimensionless constant which charac- 
terizes the noise intensity and A" is the total number of the free charge carriers or 
fluctuating particles. Since the current is proportional to the number of particles, the 
power spectral density is proportional to the number of particles A" too. Therefore, 
it is likely that every particle influences on the spectral power additively and it is 
logical to search of the 1/f fluctuations in the few-particle systems affected by some 
perturbations. 



X 



MODELS AND RESULTS 



One of the relatively universal and approved mechanisms which generates 1/f 
noise is the intermittency [2]. Intermittent chaos is often observed phenomenon in 
nonlinear dynamical systems. Investigating a transition from chaotic to nonchaotic 
behavior in the randomly driven systems [3] we have observed the intermittency 
route to chaos too [4]. Therefore, it is interesting to analyze the power spectral 
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density of the flux of small number of particles, or even one particle, moving in 
some contours and perturbed by external random forces. 

The intensity of one-dimensional flux of small particles in some space point 
x s may be expressed as 

J(t) = £sign( Ufc ) (2) 
k 

where tk is a sequence of transit times when one of the particles crosses the point 
x s and sign(f fc ) is the sign of velocity in this point. The power spectral density of 
the flux (2) is 




where T is the whole observation time and the brackets (...) denote averaging over 
the realizations. When the averaged flux I is nonzero the dominant term in the 
power spectral density at low frequencies regardless of the detailed structure of the 
flux is Sdc(f) = I /T(7rf) 2 . To eliminate this term we divide the whole observation 
time T into pieces of equal length and calculate the averaged power spectral density 
of the difference of two pieces fluxes [6] . 

We analyze the power spectral density of particles moving in some potential 
and perturbed by the random noise modelled like in papers [3, 4] by the resets 
of velocity of all particles after every time interval r according to the replacement 
v new (kT) = av old (kT) + f3v r k an ) k — 1,2, . . ..Here a and (3 are positive parameters 
and v r k an are random velocities from the Gaussian distribution of variance a — 1. 
First of all we have calculated the power spectral density of the flux of particles 
bounded in the Duffing potential V (x) = x 4 —x 2 . Under some conditions a transition 
from chaotic to nonchaotic behavior and synchronization of particles in such system 
may be observed [3-5] . The calculated spectral density of the current of particles in 
such systems usually was very close to the white noise. Only for a = 1 we observe 
the 1/f 2 noise, while for a close to 1 the noise dependence on the frequency is 
between 1/f 2 and 1/f (see Fig. 1(a)). A very small component of 1/f noise in the 
flux fluctuations of such system may be explained on the basis of Eq. (1) according to 
which the intensity of 1/f noise component is proportional to the squared averaged 

— 2 

current / . But in such model the averaged current equals zero. 

Therefore further we have analyzed the power spectral density of the current 
of particles moving in the closed contour and perturbed by the common for all 
particles noise. The simplest equations of motion for such model are of the form [5] 

v = F(x) — 7f , x — v (4) 

Here F(x) = F(x + L) is the (regular) periodic in the space force with L being the 
length of the contour and 7 is the friction coefficient. Because of the directed drift the 
averaged current in this model is nonzero. We have investigated the dependence of 
the power spectral density on the nonlinearity of equations of motion and on intensity 
of the perturbation. Moreover, we have analyzed the ensemble of particles moving 
with the friction in the periodically time-dependent force F(x + L,t + T p ) = F(x, t) 
and affected of the identical for all particles noise. We have observed the relatively 
weak sensitivity of the spectral density dependence on the frequency to the nonlinear 
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and periodically time-dependent terms in the equation for the velocity. Moreover, 
the spectral density is not sensitive to the number of the particles N moving in the 
contour too. Therefore, we show results of calculations for the simplest versions 
of the model (4), i.e. for the motion of one particle affected by the constant force 
and for the free motion with perturbations modelled by the velocity replacements 
every time interval r (see Figs. 1 (b) and (d)). For some parameters r when a ^ 1 
we observe the S(f) dependence on the frequency / (in same range of frequencies) 
close to the 1 / /-dependence. Although, when we analyze the current power spectral 
density in the large range of frequency we can, as a rule, discover the Lorentzian- 
type curves. However the increase of the S(f) with the decrease of the frequency 
usually is not like 1/ f 2 but often as 1/ /. 

From the analyzes above it appears that the current power spectral density in 
the most degree depends not on the type of the potential (and force F(x, t)) in which 
the particles move but much more on the type and intensity of the perturbation. 
Therefore, we have also analyzed the simplest model for the perturbed flux: the 
free motion of the particles with the Brownian-type additional fluctuations of the 
transition times. For such model the transit times may be calculated from the 
iterative sequence tk = tk-i + 1 + (3b~t r k an . Here the period of the free motion in 
the closed contour is chosen to be equal 1 and 5t r k an are random time increments 
from the Gaussian distribution of variance a — 1. Results of calculations of S(f) 
for different values of intensity of perturbation (5 are presented in Fig. 1 (c). With 
increase of the perturbation intensity we observe transition from the white noise at 
P < 1 to the Lorentzian-type, 1/ /- and l// 1 - 3 -type noise at f3 > 2. 

Summarizing, we have analyzed the problem of necessary and sufficient 
conditions for appearance of the 1 / /-type fluctuations in the current of the simple 
systems of particles moving in some contour. Even in such simplest models consisting 
of few or even one particle moving in the closed contours and perturbed by the 
appropriate external noise we may observe the power spectral density S(f), which 
behaves like 1/ f 5 (S ~ 1) in same range of frequencies. 
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Fig. 1 . Power spectral densities as functions of frequency for flux of particle: (a) 
in the Duffing potential with x = 3, p = 1 and different a; for motion in closed 
contour according to Eq.(4): (b) with F(x) = 1, 7 = 0.1, L = 2, p = 1; (d) with 
F(x) = y=0, L = 1 , p = 1 ; and (c) for the iterative sequence of transit times. 



